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In a bilayer consisting of an insulator (I) and a ferromagnetic metal (FM), interfacial spin orbit 
scattering leads to spin mixing of the two conducting channels of the FM, which results in an uncon¬ 
ventional anisotropic magnetoresistance (AMR). We theoretically investigate the magnetotransport 
in such bilayer structures by solving the spinor Boltzmann transport equation with generalized 
Fuchs-Sondheimer boundary condition that takes into account the effect of spin orbit scattering at 
the interface. We find that the new AMR exhibits a peculiar angular dependence which can serve as 
a genuine experimental signature. We also determine the dependence of the AMR on film thickness 
as well as spin polarization of the FM. 
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Anisotropic magnetoresistance (AMR) is a generic 
magnetotransport property of ferromagnetic metals. In 
general, the longitudinal resistance of a bulk polycrys¬ 
talline ferromagnetic metal only depends on the relative 
orientations of the magnetization vector and the current, 
which can be cast in the form [I], 

P = Po + A Pb(ie ■ m ) 2 (1) 

where j e = j e /je is the unit vector in the direction of the 
current density, m is the unit vector in the direction of 
the magnetization, p 0 is the isotropic longitudinal resis¬ 
tivity and A p b quantifies the magnitude of the bulk AMR 
effect (typically A p b ~ l%for transition metals and their 
alloys). The effect has found many practical applications 
in magnetic recording and sensor devices [2] ■ 

Recently, the AMR effect has also played a key role 
in measurements of spin Hall angle EH5] as well as 
spin torque generation EHg in FM/heavy-metal and 
FM/topological-insulator (TI) bilayers. The structural 
inversion asymmetry of these structures, combined with 
strong spin orbit coupling in the non-magnetic layers, 
generates a large Rashba-type spin-orbit coupling da¬ 
ng 

V s . 0 . = -\ 2 C V' (z) a • (p x z) (2) 

where er is the Pauli spin matrix, p is the momentum op¬ 
erator, A c is the effective Compton wave length, z is the 
unit vector normal to the interface, V(z) is the potential 
in the vicinity of the interface (which only varies in the 
direction), and V'{z) is its derivative, which is large 
only in the interfacial region. A natural question arises: 
will the interfacial spin orbit interaction alter the AMR 
in the FM layer? At first glance, one might think the 
spin orbit interaction, commuting with the total in-plane 
momentum, p x and p y , cannot alter the in-plane resis¬ 
tivity. However, this argument fails in a ferromagnetic 
metal since the in-plane momenta of either spin com¬ 
ponent are not separately conserved, and the spin-orbit 



FIG. 1: Schematics of the transverse AMR effect induced by 
surface spin orbit scattering in a FM/I bilayer. The longitu¬ 
dinal resistivity changes when the magnetization m is rotated 
in the plane perpendicular to the electric field E. Specifically, 
the resistivity is at a maximum when m is perpendicular to 
the interface (i.e., (5 = 0) and reaches a minimum when m 
lies in the plane of the interface (i.e., /3 = ±7r/2). This un¬ 
usual AMR effect arises from spin mixing of the conducting 
channels of the FM, which depends on the relative directions 
of the magnetization m and the effective spin orbital field 
H so ~zxE. 


coupling transfers momentum from one spin channel into 
the other. 

In this paper, we show that, in the presence of surface 
spin orbit scattering, the AMR of a ferromagnet exhibits 
an angular dependence that is distinctly different from 
the conventional one described by Eq. (1). As shown in 
Fig. 1, when the magnetization vector m is swept in the 
plane perpendicular to the applied electric field E, a vari¬ 
ation in the longitudinal resistivity occurs, which has no 
analogue in Eq. 0 : the resistivity has a maximum when 
m is along the z— axis (i.e., normal to the film plane) 
and reaches a minimum when m is along the y— axis 
(i.e., orthogonal to both E and z), even though the an¬ 
gle between the magnetization and the current does not 
change. The physical origin of this unconventional an- 
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gular dependence lies in the concerted actions of surface 
spin orbit scattering and spin asymmetry in the conduc¬ 
tivity of the FM, which can be understood qualitatively 
within the two-current model mm- The surface spin 
orbit scattering plays a crucial role in mixing the two 
parallel current channels of majority and minority spins; 
moreover, the degree of spin mixing depends on the rel¬ 
ative orientation of the magnetization and the effective 
magnetic field H so ~ j e x z seen by the electron spin. 
Specifically, spin mixing is strong when m is perpendic¬ 
ular to H so while it is weak when m is aligned with H so . 
The spin mixing causes, as long as the resistivities of the 
two current channels are not identical, a redistribution 
of current by decreasing the resistivity of the channel 
with higher resistivity and increasing the resistivity of 
the channel with lower resistivity: this results in an over¬ 
all increase of the total resistivity mEi. The largest 
resistivity therefore coincides with the largest degree of 
spin mixing, which occurs when m is perpendicular to 

n so . 

In the remainder of this paper we present a quantita¬ 
tive theory of the spin-orbit driven AMR in a FM/I bi¬ 
layer with strong interfacial spin orbit coupling. At vari¬ 
ance with previous studies Il5( - ll9j we exclude from our 
consideration heavy metals with strong spin-orbit cou¬ 
pling: in fact, we assume that the spin-orbit coupling is 
negligible in the bulk of the ferromagnet, so as to avoid 
contamination from the conventional bulk AMR effect. 

Our theoretical analysis is based on the spinor form 
of the semiclassical Boltzmann equation: the non- 
equilibrium distribution function, /(k.z), is a 2 x 2 ma¬ 
trix in spin space H3 HU- For simplicity, we assume 
same Fermi wave vector but different relaxation times 
for majority and minority electrons, just as in the sem¬ 
inal paper by Valet and Fert in calculating the CPP- 
GMR [22]. Such simplification is justified since the 
essence of our effect is in the difference of the relaxation 
times (and hence the resistivities) of majority and minor¬ 
ity spins, while the exchange splitting is responsible for 
spin dephasing of the transverse spin component [201128] . 

The equation of motion for /(k,z) in the steady state 
is 


df(K z ) w fdf 0 (k)\ f 1 f—i x t n G 

'■Tr- tS '-lT t rJ' = -2{ T 

( 3 ) 

where fo(k) is the equilibrium distribution function, 


5.-1 


(t 7 ) 1 (7 - per • mj is the spin dependent re¬ 
laxation rate with (r 7 ) -1 = + (tj .) -1 /2 and 


P = 


( r i) 1 - Or) 1 / ( T t) 1 + ( r l) 1 bein g> respec¬ 


tively, the average momentum relaxation rate and the 
spin asymmetry in resistivity, r-j- and tj, being the mo¬ 
mentum lifetimes of the two spin channels, and { ,} 
standing for an anticommutator. Notice that in the col¬ 
lision term of Eq. ([3]) we require that the difference 

Sf{k„z) = f{k,z) - f Le .(k,z) (4) 


between the non-equilibrium distribution and a local 
equilibrium distribution |24H26j 


fi.e.(k,z) = fo(k)i+^ 
OEk 


Mo (z) I + er ■ AG 0 ) ( 5 ) 


tend to zero for long times. The parameters p 0 (z) and 
fz s (z) of the local distribution are fixed in such a way 
that the condition 


/ 


d 3 k 


f{k,z) - fi. e .(k,z) 


= 0 


( 6 ) 


is self-consistently satisfied. By doing this, we satisfy the 
physical requirements of particle and spin conservation in 
the collision processes, as well as the continuity equations 
that go with them. 

To solve the Boltzmann equation, we further separate 
the distribution function into an equilibrium part and a 
small non-equilibrium perturbation, i.e., 


/(M = /o( fc 7+§| 


9 (M) / + h(k,z) 


( 7 ) 


where g (k,z) I and h(k,z) • <x are the spin-independent 
and spin-dependent components of the non-equilibrium 
distribution. By inserting Eq. 0 into Eq. 0 , we obtain 
a set of coupled equations for the scalar and vector parts 
of the distribution function 


d 9 p 

Vz a- eE x V x 

oz 


5-Mo - P 




( 8 ) 


and 


9h\\ h\\ -/Gy -p(s-Mo) 

Vz dz t' 


( 9 ) 


<9hj_ _ hj. - fi s± 

>z dz ~ t' 


where h\\ = m-h(k,z), = m -fJ, s (z), hj_ = 

(m x h) xm and fi s± = (m x/i s ) xm. The equations 
([8|)-([l0| have the general solutions 
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g ± (k,z) = erv x E x + A±e T 


p-p)* 

\ v z\ r ' 


± „T 


B k e 


(l+p)z 

\vz\t' 


E 


dt 


do (t) + a/x s |, (t) 


- (!-< 


J dt 6 


pH- 

•z\r> 


( 11 ) 


/ijf (k,z) = perv x E x + A±e T — 


ci-p)» 


4^ 


(i+p)* 


and 

h± (k,z)=C±e T ^ 


dtfi a± (1) , (13) 


where the superscript + labels the solution for v z > 0 
and the subscript — for v z < 0. The sum over a runs 
over the values a = ±1. The four unknown constants 
Ak, -Bki and Ck (where Ck is a vector orthogonal to m, 
hence with only two components) will now be determined 
from the boundary conditions. 

Up to this point, the interfacial spin-orbit interaction 
has not appeared in our calculations. In particular, the 
collision term in Eq. ([3]) did not contain it, and therefore 
conserved spin. The spin-orbit coupling appears in the 
boundary condition that connects the distribution func¬ 
tion for electrons impinging on the interface (label —) to 
the distribution function for electrons that are scattered 
off the interface (label +). Specifically, in order to take 
into account the rotation of spin upon scattering off the 
interface with the potential 


V scat . = V b Q (-z) - (V b X 2 c ) 8{z)a ■{p x z) (14) 


(where V b is the barrier height of the insulator, 0 (z) is 
the unit step function, z is the unit vector normal to the 
interface and the delta function confines the SO coupling 
to the interface at z = 0 ), we introduce the following 
spinor generalization of the Fuchs-Sondheimer boundary 
condition [27H25] : 


/+(k; 2 = 0) = sjftf-(k;z = 0 )R + (1 - S/ ) /„/ (15) 


where the superscripts + and — correspond to the dis¬ 
tribution functions with v z > 0 and v x < 0 respectively, 
the parameter s/ varies between 0 and 1 , characterizing 
the fraction of electrons being specularly reflected [29] 
(sj = 1 when the interface is perfectly smooth and sj = 0 
when the interface is extremely rough) and R is a 2 x 2 
reflection amplitude matrix in spin space which captures 
the spin rotation of the reflected electrons. Explicitly, 



(A c k b ) 4 q 2 I + 2i (X c k b ) z k z cr • (q x z) 


(k — ik z ) - (X c k b )'q 2 


(16) 


where k = \/k 2 — k 2 with k b = 2m*V b /h 2 . The deriva¬ 
tion of R is presented in the appendix. It must be pointed 




a 


dt 


do (t) + ads II (t) 


.(l-qp)(z-/ 


dt 6 


( 12 ) 


out that the boundary condition described by Eq. (15) 


neglects the interference between incident and reflected 
states. However, we neglect this quantum coherence by 
assuming a rough interface. This is justified when the 
characteristic size as well as the correlation length of the 
surface roughness is comparable to the Fermi wavelength 
in which case the phase coherence is destroyed by surface 
roughness 130 . 

For simplicity, we assume spin independent specular 
reflection only at the outer surface of the FM layer (z = 
d ), i-e., 


f + (k-z = d) = f (k \z = d) 


(17) 


Neglecting spin-dependent scattering from the other sur¬ 
face simplifies the calculation without altering any qual¬ 
itative features of the results. By inserting Eq. 0 into 
the boundary conditions as well as Eq. ([ 6 ]) , we can find 
unique solutions for /(k,z) and the charge current den¬ 
sity can be calculated as 

je W = 7XrA Tr f dk /( k ^) v ( 18 ) 

yZlT J J 

After some algebra, we find the charge current den¬ 
sity up to second order in the spin orbit coupling, i.e., 

O ((A c k b f) 


je (z) = C 0 E X | [l-aW ( 2 ) - 4x ( 2 ) 


x+a yx ( 2 ) y j 


(19) 

where Co = e 2 Tkp/3nml is Drude conductivity and two 
position dependent coefficients read 

4* (~) = (1 ~ si) (! + VP) f p° ( z ) ( 20 ) 


4* ( z ) = sip (A c k F ) 4 [4 - (m 2 + 3m 2 )] ^ ctG p<j ( 2 ) , 

( 21 ) 


and 


a yx ( 2 ) = Sjp (A c k F ) 4 m x m y ^ aG pa ( 2 ) , (22) 


where 


F piJ (z) = ^ j dt; 


3 f 1 (l-C 2 )cosh 


( l-crp)(d-z ) 
JA 0 (1-P 2 ) 


exp 


(1 — ap)d 

— si exp 

(1 — crp)d 

«Ao(l -p 2 ) 

?Ao(l—P 2 ) 



(23) 


































4 



FIG. 2: The transverse AMR A p[ a J as a function of the spin 
polarization p for several different s/. Parameters: XJcf = 
0.05, d/A 0 = 1. 


and 


GW 


£(1 — £ 2 ) 3 / 2 cosh 


(l-gp)(d-z) 
«Ao(l—P 2 ) 


exp 


(1 — ap)d 


(1 —<rp)d 

|_«Ao(1-P 2 )J 

s j exp 

£A 0 (l-p 2 ) 


(24) 

and Aq = vp(rf + T^)/2 is the average electron mean free 


path. The first term, 0 : 2 , is independent of the magne¬ 
tization direction: it is the resistivity due to interfacial 
roughness [3D). The third term, a yx , corresponds to the 
well-known planar Hall effect mm- 

The second term, a^J, describes the new AMR effect. 
We note that this effect is of second order in the spin- 
orbit coupling and vanishes when the spin polarization 
p = 0. The experimentally relevant quantity is the spa¬ 
tially averaged longitudinal resistivity, which is obtained 
from the formula p~ x (d) = (1/d) Jq dz j e , x (z) /E x . As 
we discussed earlier, the bulk spin-orbit coupling (not in¬ 
cluded in our calculation) produces a conventional AMR 
with angular dependence shown in Eq. (1). Therefore in 
general, the longitudinal resistivity of the FM thin film 
should take the form 


Pxx ( d ) =Po + A P { so +s)m l - A P ( so m l (25) 

where the first term on the rhs of Eq. ( [25] ) is the isotropic 
resistivity and the second term is the AMR with conven¬ 
tional angular dependence of (j e • m) 2 to which both bulk 
and surface spin orbit coupling may contribute. The most 
interesting term is the third term which is solely due to 
the surface spin orbit scattering and can be distinguished 
from the bulk AMR based on the different angular de¬ 
pendence. 

In Fig. 2, we show A pi s J (normalized by p 0 ) as a func¬ 
tion of the spin polarization p. This figure delivers two 
main messages. First, A plo is positive when p / 0. 
This confirms the angular dependence of the AMR we 
sketched in Fig. 1(b). The second message is that A p[ s J 



FIG. 3: The transverse AMR Ap// as a function of d/Xo for 
several different p. Parameters: X c k,F = 0.05 and sj = 0.8. 


is an even function of the spin polarization p. This is 
consistent with our two-current model: the spin mixing 
resistivity only relies on the absolute value of the differ¬ 
ence between the two conduction channels but not on the 
sign of that difference. 


In Fig. 3, we show the thickness dependence of the 
A pso for several values of spin polarization. When the 
FM layer thickness is much larger than the mean free 
path, i.e., d Ao, A p{ s J exhibits a standard 1/d thick¬ 
ness dependence as can be analytically worked out via 
Eq. m. 


Lastly, let us consider the choice of materials for the 
observation of our predicted AMR. In order to obtain a 
sizable transverse AMR, it is essential to have a FM/I bi¬ 
layer with a large difference between the conductivities of 
majority and minority spin carriers in the ferromagnetic 
metal, and a strong spin-orbit interaction at the FM/I 
interface. A very promising system in this respect is Py 
grown on top of a TI such as Bi 2 Se 3 , Bi 1 . 5 Sbo. 5 Te 1 . 7 Se 1.3 
or Sn-doped Bi 2 Tei. 7 Sei . 3 [311 - 153] . Recently, large spin 
transfer torque |S] and spin-charge conversion [M| effects 
were observed in these FM/TI bilayers, indicating the 
presence of strong spin orbit coupling at the interface. 
Non-topological-oxide/ferromagnetic interfaces may also 
provide large spin-orbit interaction, as implied by tun¬ 
neling AMR studies in Fe/MgO/Fe junctions [35] [SB] 
and by magnetic anisotropy analysis at AlOx/Co inter¬ 
face E3I33. As a final point, we note that the surface 
spin orbit scattering mechanism should, in principle, also 
contribute to the transverse AMR that was previously 
found in Pt/Co/Pt [38] and Py/YIG [3D] layered sys¬ 
tems. 
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Appendix: Spinor form reflection amplitude in the Now we are ready to determine R and T matrices by 

presence of interface Rashba spin orbit coupling the following quantum mechanical matching conditions 


In this appendix, we derive the spinor form reflection 
amplitude given by Eq. (16) in the main text. First, we 
write down the following piece-wise scattering wave func¬ 
tions corresponding to the interfacial potential described 
by Eq. (14) in the main text 

0k (r\\,z > 0) = _L e -^ V q ' r i, x+ _Le ifc ‘V q ' r iiI?x 

(Al) 

and 

V> k (r,;z<0) = ^e K V qr lif X (A2) 


0k ( r lb 0+ ) = V>k ( r |b0~ 


(A3) 


and 


0k ( r |b 0 + )—0k ( r ||,0 )= (fc b A c ) 2 er ■ (p x z) 0 k (r||,O) 

(Al) 

By placing the scattering wave functions into the above 
two equations, we find 


(-T + R) X = T X 


(A5) 


and 


where ry is the in-plane position vector, q and k z are the 
in-plane and perpendicular-to-plane wave vectors respec¬ 
tively, k = \Jkl — fc?with kb = \J2m*Vb/h 2 , R and T are 
the 2 x 2 reflection and transmission amplitude matrices 
in spin space, and x is an arbitrary spinor. 


— kT — ik z I + ik z R] x = 


(kb A c ) 2 (r ■ (q x z) 


T X 

(A6) 


Combining the two equations, we find an equation for R 
only, 


| [(k - ik z ) / + (k b A c ) 2 <t • (q x z) R + ( ik z + k) I + (k b A c ) 2 er • (q x z) j X = 0 
For any X , the equation is valid if 

(k — ik z ) i + (k b A c ) 2 <j • (q x z) R + ( ik z + k) I + (k b A c ) 2 <r • (q x z) = 0 


(A7) 


(AS) 


From Eq. (A 8 ), we find the reflection amplitude matrix 


R = 


-k b + (fc b A c ) 4 q 2 I + 2 i (k b \ c Y k z cr ■ (q x z) 


(ik z — k) - (k b X c )' q 2 


(A9) 

. Note that R is a unitary matrix as can be easily check 
via Eq. (A91. This unitarity is required by flux conser¬ 
vation. 


* Electronic address: zhangshule@missouri.edu 
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